Relaying method for protecting transformers

ABSTRACT

The present invention relates to a relaying method using the ratio of induced voltages or the ratio of flux linkage increments. The protective relaying method for power transformers with one or more phases includes the first step of obtaining primary and secondary voltages and currents of the transformer; the second step of calculating induced voltages, induced voltage differences, ratio of primary and secondary induced voltages, or ratio of primary and secondary induced voltage differences from the currents and the voltages; the third step of calculating at least one predetermined decision parameter derived from at least one predetermined equation; the fourth step of deciding whether an internal winding fault occurs by comparing the decision parameter to the induced voltages, the induced voltage differences, the ratio of primary and secondary induced voltages, or the ratio of primary and secondary induced voltage differences.

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates to a relaying method for transformer protection. More specifically, the present invention relates to a relaying method using the ratio of induced voltages or the ratio of flux linkage increments.

2. Description of the Related Art

Current differential relays with harmonic constraints are widely used for power transformers. These relays determines a fault if a differential current between the primary and secondary windings is larger than a threshold value. The relays may malfunction in case of magnetic inrush or overexcitation. In order to prevent malfunction they use the second harmonic, third harmonic, or the total harmonic component of the current.

However, because of increase of operating voltages and length of transmission line, the differential currents during internal transformer faults can contain large harmonic components. Therefore, the security of the differential relays using harmonic restraints is a matter of concern. Also, these relays are inherently insensitive to turn-to-turn faults.

Recently, some algorithms use electromagnetic equations of transformers.

One of them is based on an equivalent circuit composed of inverse inductance. Another one is a flux-restrained current-differential relay. Both of them use winding currents those are practically unavailable in some situations such as a three-phase wye-delta transformer. Some elements used in the first method are nonlinear depending on the operating condition and thus it is very difficult to determine their values. In case of the second method, the winding resistance is neglected. But, practically it cannot be neglected. Thus, if there exits a winding resistance, large errors are contained.

SUMMARY OF THE INVENTION

Thus, the invention intends to overcome the above-mentioned problems and the purpose of the invention is to provide a protective relaying method for power transformers. In other words, the ratio of primary and secondary flux linkage increments is the same as the turn ratio in case of magnetic inrush and overexcitation, while the ratio is different from the turn ratio in case of internal winding faults. Thus, this invention provides a protective relaying method for power transformers; the method calculates the ratio of primary and secondary flux linkage increments from the primary and secondary voltages and currents and then compares the ratio with the turn ratio to discriminate magnetic inrush and internal winding faults.

And, another purpose of the invention is to provide a protective relaying method for power transformers using the ratio of primary and secondary induced voltages. In other words, the ratio of induced voltages is the same as the turn ratio in case of magnetic inrush and overexcitation, while the ratio is different from the turn ratio in case of internal winding faults. Thus, this invention provides a protective relaying method for power transformers; the method calculates the ratio of primary and secondary induced voltages from the primary and secondary voltages and currents and then compares the ratio with the turn ratio to discriminate magnetic inrush and internal winding faults.

According to a feature of the present invention to achieve the above object, the protective relaying method for power transformers with one or more phases includes the first step of obtaining primary and secondary voltages and currents of said transformer; the second step of calculating induced voltages, induced voltage differences, ratio of primary and secondary induced voltages, or ratio of primary and secondary induced voltage differences from said currents and said voltages; the third step of calculating at least one predetermined decision parameter derived from at least one predetermined equation; the fourth step of deciding whether an internal winding fault occurs by comparing said decision parameter to said induced voltages, said induced voltage differences, said ratio of primary and secondary induced voltages, or said ratio of primary and secondary induced voltage differences.

According to another feature of the present invention to achieve the above object, the protective relaying method for power transformers with one or more phases includes the first step of obtaining primary and secondary voltages and currents of said transformer; the second step of calculating flux linkage increments, flux linkage difference increments, the ratio of primary and secondary flux linkage increments, or ratio of primary and secondary flux linkage difference increments from said currents and said voltages; the third step of calculating at least one predetermined decision parameter derived from at least one predetermined equation; and the fourth step of deciding whether an internal winding fault occurs by comparing said decision parameter to said flux linkage increments, said flux linkage difference increments, said the ratio of primary and secondary flux linkage increments, or said ratio of primary and secondary flux linkage difference increments.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows the drawing of a two winding single-phase transformer.

FIG. 2 shows the drawing of a two winding three-phase Y—Y transformer.

FIG. 3 shows the drawing of a two winding three-phase Y-Δ transformer.

FIG. 4 shows the drawing of the model system including a two winding three-phase Y-Δ transformer in which the protective relaying method of this invention is applied.

FIG. 5 shows the drawing of currents of each phase in case of no remanent flux in the transformer.

FIG. 6 shows the drawing of calculated induced voltages' differences of Detector 1 among three Detectors.

FIG. 7A to FIG. 7D show the drawings of three Detectors and a trip signal in case of no remanent flux according to the first preferred embodiment of this invention.

FIG. 8A to FIG. 8C show the drawings of currents of each phase in case of +80% remanent flux of the saturation point in the transformer.

FIG. 9A to FIG. 9D show the drawings of three Detectors and a trip signal in case of +80% remanent flux according to the first preferred embodiment of this invention.

FIG. 10A to FIG. 10D show the drawings of three Detectors and a trip signal in case of the B-phase internal turn-to-earth fault at the 40% point from the neutral according to the first preferred embodiment of this invention.

FIG. 11A to FIG. 11C show the drawing of calculated flux linkage difference increments of Detector 1 among three Detectors according to preferred embodiment 2.

FIG. 12A to FIG. 12D show the drawings of three Detectors and a trip signal in case of no remanent flux according to the second preferred embodiment of this invention.

FIG. 13A to FIG. 13C show the drawings of currents of each phase in case of +80% remanent flux of the saturation point in the transformer according to preferred embodiment 2.

FIG. 14A to FIG. 14D show the drawings of three Detectors and a trip signal in case of +80% remanent flux according to the second preferred embodiment of this invention.

FIG. 15A to FIG. 15D show the drawings of three Detectors and a trip signal in case of the B-phase internal turn-to-earth fault at the 40% point from the neutral according to preferred embodiment 1 of this invention.

FIG. 16 shows the drawing of a two winding three-phase Δ-Δ transformer.

FIG. 17 shows the drawing of a two winding three-phase Δ-Y transformer.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

Preferred embodiments of the present invention will now be described in detail with reference to the accompanying drawings.

The first preferred embodiment of the invention relates to a protective relaying algorithm for power transformers using the ratio of primary and secondary induced voltages. FIG. 1 shows a two winding single-phase transformer. If there is no internal fault of the two-winding single-phase transformer, the primary and secondary voltages at any instant can be given by the following Equation 1 and Equation 2. $\begin{matrix} {v_{1} = {{R_{1}i_{1}} + {L_{l1}\frac{i_{1}}{t}} + e_{1}}} & \left\lbrack {{Equation}\quad 1} \right\rbrack \\ {v_{2} = {{{- R_{2}}i_{2}} - {L_{l2}\frac{i_{2}}{t}} + e_{2}}} & \left\lbrack {{Equation}\quad 2} \right\rbrack \end{matrix}$

where, R₁, R₂ are the primary and secondary winding resistances, respectively, and L₁₁, L₁₂ are the primary and secondary leakage inductances, respectively, and e₁, e₂ are primary, secondary induced voltages, respectively.

Rearranging Equation 1 and Equation 2 yields Equation 3 and Equation 4. $\begin{matrix} {e_{1} = {v_{1} - {R_{1}i_{1}} - {L_{l1}\frac{i_{1}}{t}}}} & \left\lbrack {{Equation}\quad 3} \right\rbrack \\ {e_{2} = {v_{2} + {R_{2}i_{2}} + {L_{l2}\frac{i_{2}}{t}}}} & \left\lbrack {{Equation}\quad 4} \right\rbrack \end{matrix}$

Thus, induced voltages e₁, e₂ can be estimated with voltages, currents, winding resistances and leakage inductances. The ratio of induced voltages (RIV) is defined by Equation 5. $\begin{matrix} {{RIV} = {\frac{e_{1}}{e_{2}} = \frac{N_{1}}{N_{2}}}} & \left\lbrack {{Equation}\quad 5} \right\rbrack \end{matrix}$

As e₁, e₂ can be estimated using Equation 3 and Equation 4 the RIV can be calculated at any instant using Equation 5. In the steady state, if there is no internal winding fault, the RIV at any instant is equal to the turn ratio, i.e. N₁/N₂ except for e₁=0 or e₂=0.

In case of magnetic inrush, as the core repeats saturation and unsaturation, a magnetizing current also repeats large and small. This phenomenon lasts for some time. However, in this case since it is not an internal winding fault, the RIV is also equal to the turn ratio even when e₁, e₂ are not sinusoidal and distorted.

On the other hand, the RIV is not the same as the turn ratio in case of internal winding faults. Thus, from this point, the first preferred embodiment of this invention detects an internal fault of a transformer based on the following fact; it is not an internal fault if the RIV is the same as the turn ratio, while it is an internal fault otherwise.

However, e₁ and e₂ are instantaneous values and the RIV can grow up to a large value when e₂ has values near zero. In the first preferred embodiment of this invention, a Detector of Equation 6 is used to detect a fault instead of Equation 5. That is, if Equation 6 is less than the threshold, it is normal; otherwise, it is an internal winding fault. $\begin{matrix} {{Detector} = {\frac{e_{1} - {\frac{N_{1}}{N_{2}}e_{2}}}{V_{1\quad {rm}\quad s}} \times 100\quad (\%)}} & \left\lbrack {{Equation}\quad 6} \right\rbrack \end{matrix}$

where, V_(1rms) is the primary rated voltage.

FIG. 2 shows the drawing of a two winding three-phase Y—Y transformer. In FIG. 2, v_(A), v_(B), v_(C) represent primary voltages and v_(a), v_(b), v_(c) means secondary voltages and i_(A), i_(B), i_(C) are primary currents and i_(a), i_(b), i_(c) are secondary currents.

In FIG. 2, voltages of the primary windings can be represented as Equation 7, Equation 8, and Equation 9. $\begin{matrix} {v_{A} = {{R_{A}i_{A}} + {L_{lA}\frac{i_{A}}{t}} + e_{A}}} & \left\lbrack {{Equation}\quad 7} \right\rbrack \\ {v_{B} = {{R_{B}i_{B}} + {L_{lB}\frac{i_{B}}{t}} + e_{B}}} & \left\lbrack {{Equation}\quad 8} \right\rbrack \\ {v_{C} = {{R_{C}i_{C}} + {L_{lC}\frac{i_{C}}{t}} + e_{C}}} & \left\lbrack {{Equation}\quad 9} \right\rbrack \end{matrix}$

where, R_(A), R_(B), R_(C) mean primary winding resistances and L_(lA), L_(lB), L_(lC) mean primary leakage inductances, and e_(A), e_(B), e_(C) mean primary induced voltages.

And, the voltages of the secondary windings are given by Equation 10, Equation 11, and Equation 12. $\begin{matrix} {v_{a} = {{{- R_{a}}i_{a}} - {L_{la}\frac{i_{a}}{t}} + e_{a}}} & \left\lbrack {{Equation}\quad 10} \right\rbrack \\ {v_{b} = {{{- R_{b}}i_{b}} - {L_{l\quad b}\frac{i_{b}}{t}} + e_{b}}} & \left\lbrack {{Equation}\quad 11} \right\rbrack \\ {v_{c} = {{{- R_{c}}i_{c}} - {L_{lc}\frac{i_{c}}{t}} + e_{c}}} & \left\lbrack {{Equation}\quad 12} \right\rbrack \end{matrix}$

where, R_(a), R_(b), R_(c) are secondary winding resistances, and L_(la), L_(lb), L_(lc) are secondary leakage inductances, and e_(a), e_(b), e_(c) are secondary induced voltages.

In case of a three-phase Y—Y transformer, all the induced voltages of the primary and secondary windings can be calculated directly from Equation 7 to Equation 12. Thus, in exactly the same manner of a single-phase transformer, three Detectors are defined by Equation 13, Equation 14, and Equation 15. $\begin{matrix} {{{Detector}\quad A} = {\frac{e_{A} - {\frac{N_{1}}{N_{2}}e_{a}}}{V_{Arms}} \times 100\quad (\%)}} & \left\lbrack {{Equation}\quad 13} \right\rbrack \\ {{{Detector}\quad B} = {\frac{e_{B} - {\frac{N_{1}}{N_{2}}e_{b}}}{V_{Brms}} \times 100\quad (\%)}} & \left\lbrack {{Equation}\quad 14} \right\rbrack \\ {{{Detector}\quad C} = {\frac{e_{C} - {\frac{N_{1}}{N_{2}}e_{c}}}{V_{Crms}} \times 100\quad (\%)}} & \left\lbrack {{Equation}\quad 15} \right\rbrack \end{matrix}$

If values of Detector A, Detector B, and Detector C are larger than the threshold, it is determined that an internal fault occurred on phase A, phase B and phase C, respectively.

The next topic is about a two-winding three-phase Y-Δ transformers.

FIG. 3 shows the connections of the primary and secondary windings of a three-phase Y-Δ transformer. In FIG. 3, v_(ab), V_(bc), v_(ca) are secondary voltages and i_(ab), i_(bc), i_(ca) are secondary currents, and i_(a), i_(b), i_(c) are line currents.

In FIG. 3, voltages of the wye connected windings are represented in Equation 7, Equation 8 and Equation 9. Equation 16, Equation 17 and Equation 18 represent voltages of delta connected windings. $\begin{matrix} {v_{ab} = {{R_{ab}i_{ab}} + {L_{lab}\frac{i_{ab}}{t}} + e_{ab}}} & \left\lbrack {{Equation}\quad 16} \right\rbrack \end{matrix}$

$\begin{matrix} {v_{bc} = {{R_{bc}i_{bc}} + {L_{lbc}\frac{i_{bc}}{t}} + e_{bc}}} & \left\lbrack {{Equation}\quad 17} \right\rbrack \\ {v_{ca} = {{R_{ca}i_{ca}} + {L_{lca}\frac{i_{ca}}{t}} + e_{ca}}} & \left\lbrack {{Equation}\quad 18} \right\rbrack \end{matrix}$

where, R_(ab), R_(bc), R_(ca) are secondary winding resistances, and L_(lab), L_(lbc), L_(lca) are secondary leakage inductances, and e_(ab), e_(bc), e_(ca) are secondary induced voltages.

As i_(ab), i_(bc), i_(ca) are not available in case of Y-Δ connection, e_(ab), e_(bc), e_(ca) cannot be calculated directly from Equation 16, Equation 17 and Equation 18. However, i_(a), i_(b), i_(c) are available. Thus, in this invention, in order to use i_(a), i_(b), i_(c) the following relations i_(ca)−i_(ab)=i_(a), i_(ab)−i_(bc)=i_(b), i_(bc)−i_(ca)=i_(c) are employed. Moreover, if there is no internal fault, R_(ab)≈R_(bc)≈R_(ca)=R and L_(lab)≈L_(lbc)≈L_(lca)=L_(l). Therefore, manipulating Equation 16, Equation 17 and Equation 18 gives the following Equation 19, Equation 20 and Equation 21. $\begin{matrix} {{e_{ca} - e_{ab}} = {v_{ca} - v_{ab} - {Ri}_{a} - {L_{l}\frac{i_{a}}{t}}}} & \left\lbrack {{Equation}\quad 19} \right\rbrack \\ {{e_{ab} - e_{bc}} = {v_{ab} - v_{bc} - {Ri}_{b} - {L_{l}\frac{i_{b}}{t}}}} & \left\lbrack {{Equation}\quad 20} \right\rbrack \\ {{e_{bc} - e_{ca}} = {v_{bc} - v_{ca} - {Ri}_{c} - {L_{l}\frac{i_{c}}{t}}}} & \left\lbrack {{Equation}\quad 21} \right\rbrack \end{matrix}$

Therefore, while e_(ab), e_(bc), e_(ca) cannot be calculated, left-hand sides of Equation 19, Equation 20 and Equation 21 can be calculated. The following Equation 22, Equation 23, and Equation 24 for the primary winding corresponding to Equation 19, Equation 20 and Equation 21 can be derived by manipulating Equation 7, Equation 8 and Equation 9. $\begin{matrix} {{e_{C} - e_{A}} = {v_{C} - v_{A} - \left( {{R_{C}i_{C}} - {R_{A}i_{A}}} \right) - \left( {{L_{lC}\frac{i_{C}}{t}} - {L_{lA}\frac{i_{A}}{t}}} \right)}} & \left\lbrack {{Equation}\quad 22} \right\rbrack \\ {{e_{A} - e_{B}} = {v_{A} - v_{B} - \left( {{R_{A}i_{A}} - {R_{B}i_{B}}} \right) - \left( {{L_{lA}\frac{i_{A}}{t}} - {L_{lB}\frac{i_{B}}{t}}} \right)}} & \left\lbrack {{Equation}\quad 23} \right\rbrack \\ {{e_{B} - e_{C}} = {v_{B} - v_{C} - \left( {{R_{B}i_{B}} - {R_{C}i_{C}}} \right) - \left( {{L_{lB}\frac{i_{B}}{t}} - {L_{lC}\frac{i_{C}}{t}}} \right)}} & \left\lbrack {{Equation}\quad 24} \right\rbrack \end{matrix}$

If there is no internal fault in case of a Y-Δ transformer, the following relationships in Equation 25 are valid. $\begin{matrix} {{\frac{e_{A}}{e_{ab}} = \frac{N_{1}}{N_{2}}},{\frac{e_{B}}{e_{bc}} = \frac{N_{1}}{N_{2}}},{\frac{e_{C}}{e_{ca}} = \frac{N_{1}}{N_{2}}}} & \left\lbrack {{Equation}\quad 25} \right\rbrack \end{matrix}$

While e_(A), e_(B), e_(C) can be calculated, e_(ab), e_(bc), e_(ca) cannot. Thus, Equation 25 cannot be applied directly to Detectors like a Y—Y transformer. Equation 26 is used for fault detection in case of a Y-Δ transformer. $\begin{matrix} {{\frac{e_{C} - e_{A}}{e_{ca} - e_{ab}} = \frac{N_{1}}{N_{2}}},{\frac{e_{A} - e_{B}}{e_{ab} - e_{bc}} = \frac{N_{1}}{N_{2}}},{\frac{e_{B} - e_{C}}{e_{bc} - e_{ca}} = \frac{N_{1}}{N_{2}}}} & \left\lbrack {{Equation}\quad 26} \right\rbrack \end{matrix}$

Equation 25 is a necessary and sufficient condition of Equation 26. Sufficient condition is proved trivially in substituting Equation 25 into Equation 26. Thus, only necessary conditions will be proved. Let three RIVs as $\begin{matrix} {{{Detector}\quad 1} = {\frac{e_{C} - e_{A} - {\frac{N_{1}}{N_{2}}\left( {e_{ca} - e_{ab}} \right)}}{V_{CArms}} \times 100\quad (\%)}} & \left\lbrack {{Equation}\quad 27} \right\rbrack \\ {{{Detector}\quad 2} = {\frac{e_{A} - e_{B} - {\frac{N_{1}}{N_{2}}\left( {e_{ab} - e_{bc}} \right)}}{V_{ABrms}} \times 100\quad (\%)}} & \left\lbrack {{Equation}\quad 28} \right\rbrack \\ {{{Detector}\quad 3} = {\frac{e_{B} - e_{C} - {\frac{N_{1}}{N_{2}}\left( {e_{bc} - e_{ca}} \right)}}{V_{BCrms}} \times 100\quad (\%)}} & \left\lbrack {{Equation}\quad 29} \right\rbrack \end{matrix}$

and ${\frac{e_{A}}{e_{ab}} = \alpha},{\frac{e_{B}}{e_{bc}} = \beta}$

Substituting the three RIVs into Equation 26 yields the following three equations $\frac{e_{C}}{e_{ca}} = {\gamma.}$

and ${{{\gamma \quad e_{ca}} - {\alpha \quad e_{ab}}} = {\frac{N_{1}}{N_{2}}\left( {e_{ca} - e_{ab}} \right)}},{{{\alpha \quad e_{ab}} - {\beta \quad e_{bc}}} = {\frac{N_{1}}{N_{2}}\left( {e_{ab} - e_{bc}} \right)}}$

In order for the former three Equations to be valid for all e_(ab), e_(bc), e_(ca), the condition of α=β=γ=N₁/N₂ should be satisfied. Thus, ${{\beta \quad e_{bc}} - {\gamma \quad e_{ca}}} = {\frac{N_{1}}{N_{2}}{\left( {e_{bc} - e_{ca}} \right).}}$

Therefore, Equation 25 is a necessary and sufficient condition of Equation 26.

Thus, in the first preferred embodiment of the invention, the Detectors for a Y-Δ transformer are given by Equation 27, Equation 28, and Equation 29. $\frac{e_{A}}{e_{ab}} = {\frac{e_{B}}{e_{bc}} = {\frac{e_{C}}{e_{ca}} = {\frac{N_{1}}{N_{2}}.}}}$

Table 1 shows the rule of detecting internal faults and the faulted phase from the Detectors of Equation 27, Equation 28, and Equation 29 in case of a Y-Δ transformer.

TABLE 1 Detector 1 Detector 2 Detector 3 Fault detection < < < No internal fault Th. Th. Th. Internal winding fault ≧ ≧ < A Th. Th. Th. phase < ≧ ≧ B Th. Th. Th. phase ≧ < ≧ C Th. Th. Th. phase ≧ ≧ ≧ Internal winding Th. Th. Th. fault

If all the three Detectors are less than threshold (Th.), it is not an internal fault. If one Detector is less than threshold and the other two Detectors are not, the faulted phase can be identified as shown in Table 1. If two or more phases are faulted, fault detection is possible. However, the faulted phase cannot be identified.

Case studies of the first preferred embodiment of a protective relaying method are now provided.

FIG. 4 shows the model system studied in this invention. A two winding three-phase Y-Δ transformer (154 kV/14.7 kV, 100 (MVA)) is used to generate fault data and inrush data using EMTP. An internal winding fault modeling method is used to model the turn-to-ground and turn-to-turn winding faults. The sampling rate of 16 sample/cycle is employed.

First case is magnetic inrush.

We will explain the two magnetic inrush cases of with 0(%) and +80(%) remanent fluxes of the saturation point in case of 0 (deg) energization angle and no load.

FIG. 5 shows three-phase currents in case of magnetic inrush with a 0(%) remanent flux, 0 (deg) energization angle, and no load.

FIG. 6 shows calculated induced voltages' differences of Detector 1 i.e. e_(C)−e_(A) and N₁/N₂(e_(ca)−e_(ab)) among three Detectors. The result indicates that the two calculated induced voltage differences are nearly the same even if inrush currents contain harmonic components.

FIG. 7A to FIG. 7D show the drawings of three Detectors and a trip signal in case of no remanent flux according to preferred embodiment 1 of this invention. Trip signal in FIG. 7, means the signal to initiate the circuit breaker after detection of a fault.

In this preferred embodiment, if a Detector is greater than 5(%) a counter is increased by one and otherwise decreased by one. In addition to that, if the counter is less than zero, it is reset to zero. If a counter exceeds four i.e. ¼ cycle, the final trip is issued. As the three Detectors are less than 5(%) the trip signal is blocked even the A-phase inrush current increases up to nearly 900 (A) after the energization.

The errors shown in Detectors of FIG. 7 are caused from approximation of differentiation terms in Equation 19 to Equation 24. In this preferred embodiment, the sampling rate of 16 (c/s) is used. However, these errors can be decreased if a higher sampling rate is used.

FIG. 8 and FIG. 9 show currents of each phase, the three Detectors, and a trip signal in case of +80% remanent flux of the saturation point in the transformer. As expected, due to the maximum flux, large inrush current (nearly 1600 (A)) is drawn. However, as shown in FIG. 9, all the three Detectors are less than 5(%) even when very large inrush currents and the trip signal is blocked.

An internal winding fault is now explained.

FIG. 10 shows the three Detectors and a trip signal in case of the B-phase internal turn-to-earth fault at the 40% point from the neutral according to preferred embodiment 1 of this invention. As expected, in case of an internal winding fault, the ratio of induced voltage differences is not the same as the turn ratio. As shown in FIG. 10, the proposed method detects the internal winding fault 4.167 (ms) after a fault.

Finally, preferred embodiment 1 proposes a transformer protective relaying algorithm using the ratio of induced voltages. The proposed algorithm calculates the induced voltages of the primary and secondary windings from the currents and voltages and then clearly discriminates between internal winding faults and magnetic inrush with the ratio. Various test results clearly demonstrate that the algorithm does not malfunction even in case of inrush with the maximum remanent flux. Moreover, it detects internal winding faults in only ¼ cycle after a fault.

Although the invention is described herein with reference to the preferred embodiment, it should be clearly understood that many variations and/or modifications of the basic inventive concepts herein taught which may be apparent to those skilled in the present art will remain within the spirit and scope of the present invention, as defined in the appended claims.

The cases of three-phase Δ-Y transformer and Δ-Δ transformer are briefly described below as examples of the above mentioned modifications.

FIG. 17 shows a two-winding three-phase Δ-Y transformer. Voltages of the primary windings can be represented by the following Equations. $\begin{matrix} {v_{AB} = {{R_{AB}i_{AB}} + {L_{lAB}\frac{i_{AB}}{t}} + e_{AB}}} & \left\lbrack {{Equation}\quad 30} \right\rbrack \\ {v_{BC} = {{R_{BC}i_{BC}} + {L_{lBC}\frac{i_{BC}}{t}} + e_{BC}}} & \left\lbrack {{Equation}\quad 31} \right\rbrack \\ {v_{CA} = {{R_{CA}i_{CA}} + {L_{lCA}\frac{i_{CA}}{t}} + e_{CA}}} & \left\lbrack {{Equation}\quad 32} \right\rbrack \end{matrix}$

Also, voltages of the secondary windings can be represented by the following Equations. $\begin{matrix} {v_{a} = {{{- R_{a}}i_{a}} - {L_{la}\frac{i_{a}}{t}} + e_{a}}} & \left\lbrack {{Equation}\quad 33} \right\rbrack \\ {v_{b} = {{{- R_{b}}i_{b}} - {L_{l\quad b}\frac{i_{b}}{t}} + e_{b}}} & \left\lbrack {{Equation}\quad 34} \right\rbrack \\ {v_{c} = {{{- R_{c}}i_{c}} - {L_{lc}\frac{i_{c}}{t}} + e_{c}}} & \left\lbrack {{Equation}\quad 35} \right\rbrack \end{matrix}$

In a similar manner to a Y-Δ transformer, the relationships of r_(AB)≈r_(BC)≈r_(CA)≈r, L_(lAB)≈L_(lBC)≈L_(lCA)≈L_(l), i_(CA)−i_(AB)=−i_(A), i_(AB)−i_(BC)=−i_(B), i_(BC)−i_(CA)=−i_(C) are used. From the above Equations, the differences of induced voltages of the primary windings can be obtained as follows: $\begin{matrix} {{e_{CA} - e_{AB}} = {v_{CA} - v_{AB} + {Ri}_{A} + {L_{l}\frac{i_{A}}{t}}}} & \left\lbrack {{Equation}\quad 36} \right\rbrack \\ {{e_{AB} - e_{BC}} = {v_{AB} - v_{BC} + {Ri}_{B} + {L_{l}\frac{i_{B}}{t}}}} & \left\lbrack {{Equation}\quad 37} \right\rbrack \\ {{e_{BC} - e_{CA}} = {v_{BC} - v_{CA} + {Ri}_{C} + {L_{l}\frac{i_{C}}{t}}}} & \left\lbrack {{Equation}\quad 38} \right\rbrack \end{matrix}$

Also, the differences of induced voltages of the secondary windings can be given by the following Equations. $\begin{matrix} {{e_{c} - e_{a}} = {v_{c} - v_{a} + \left( {{R_{c}i_{c}} - {R_{a}i_{a}}} \right) + \left( {{L_{lc}\frac{i_{c}}{t}} - {L_{la}\frac{i_{a}}{t}}} \right)}} & \left\lbrack {{Equation}\quad 39} \right\rbrack \\ {{e_{a} - e_{b}} = {v_{a} - v_{b} + \left( {{R_{a}i_{a}} - {R_{b}i_{b}}} \right) + \left( {{L_{la}\frac{i_{a}}{t}} - {L_{l\quad b}\frac{i_{b}}{t}}} \right)}} & \left\lbrack {{Equation}\quad 40} \right\rbrack \\ {{e_{b} - e_{c}} = {v_{b} - v_{c} + \left( {{R_{b}i_{b}} - {R_{c}i_{c}}} \right) + \left( {{L_{l\quad b}\frac{i_{b}}{t}} - {L_{lc}\frac{i_{c}}{t}}} \right)}} & \left\lbrack {{Equation}\quad 41} \right\rbrack \end{matrix}$

Like a Y-Δ transformer, Equation 42 is a necessary and sufficient condition of Equation 43. From the above fact, the three Detectors of Equation 44, Equation 45, and Equation 46 can be used as the predetermined decision parameter. $\begin{matrix} {{\frac{e_{CA}}{e_{c}} = \frac{N_{1}}{N_{2}}},{\frac{e_{AB}}{e_{a}} = \frac{N_{1}}{N_{2}}},{\frac{e_{BC}}{e_{b}} = \frac{N_{1}}{N_{2}}}} & \left\lbrack {{Equation}\quad 42} \right\rbrack \\ {{\frac{e_{CA} - e_{AB}}{e_{c} - e_{a}} = \frac{N_{1}}{N_{2}}},{\frac{e_{AB} - e_{BC}}{e_{a} - e_{b}} = \frac{N_{1}}{N_{2}}},{\frac{e_{BC} - e_{CA}}{e_{b} - e_{c}} = \frac{N_{1}}{N_{2}}}} & \left\lbrack {{Equation}\quad 43} \right\rbrack \\ {{{Detector}\quad 1} = {\frac{e_{CA} - e_{AB} - {\frac{N_{1}}{N_{2}}\left( {e_{c} - e_{a}} \right)}}{V_{CArms}} \times 100\quad (\%)}} & \left\lbrack {{Equation}\quad 44} \right\rbrack \\ {{{Detector}\quad 2} = {\frac{e_{AB} - e_{BC} - {\frac{N_{1}}{N_{2}}\left( {e_{a} - e_{b}} \right)}}{V_{ABrms}} \times 100\quad (\%)}} & \left\lbrack {{Equation}\quad 45} \right\rbrack \\ {{{Detector}\quad 3} = {\frac{e_{BC} - e_{CA} - {\frac{N_{1}}{N_{2}}\left( {e_{b} - e_{c}} \right)}}{V_{BCrms}} \times 100\quad (\%)}} & \left\lbrack {{Equation}\quad 46} \right\rbrack \end{matrix}$

FIG. 16 shows a two-winding three-phase Δ-Δ transformer. Voltages of the primary windings can be represented by the following Equations. $\begin{matrix} {v_{AB} = {{r_{AB}i_{AB}} + {L_{lAB}\frac{i_{AB}}{t}} + e_{AB}}} & \left\lbrack {{Equation}\quad 47} \right\rbrack \\ {v_{BC} = {{r_{BC}i_{BC}} + {L_{lBC}\frac{i_{BC}}{t}} + e_{BC}}} & \left\lbrack {{Equation}\quad 48} \right\rbrack \\ {v_{CA} = {{r_{CA}i_{CA}} + {L_{lCA}\frac{i_{CA}}{t}} + e_{CA}}} & \left\lbrack {{Equation}\quad 49} \right\rbrack \end{matrix}$

Also, voltages of the secondary windings can be represented by the following Equations. $\begin{matrix} {v_{ab} = {{r_{ab}i_{ab}} + {L_{lab}\frac{i_{ab}}{t}} + e_{ab}}} & \left\lbrack {{Equation}\quad 50} \right\rbrack \\ {v_{bc} = {{r_{bc}i_{bc}} + {L_{lbc}\frac{i_{bc}}{t}} + e_{bc}}} & \left\lbrack {{Equation}\quad 51} \right\rbrack \\ {v_{ca} = {{r_{ca}i_{ca}} + {L_{lca}\frac{i_{ca}}{t}} + e_{ca}}} & \left\lbrack {{Equation}\quad 52} \right\rbrack \end{matrix}$

In a similar manner to a Δ-Y transformer, the relationships of r_(AB)≈r_(BC)≈r_(CA)≈r₁, L_(lAB)≈L_(lBC)≈L_(lCA)≈L_(l1), i_(CA)−i_(AB)=−i_(A), i_(AB)−i_(BC)=−i_(B), i_(BC)−i_(CA)=−i_(C), r_(ab)≈r_(bc)≈r_(ca)≈r₂, L_(lab)≈L_(lbc)≈L_(l2), i_(ca)−i_(ab)=i_(a), i_(ab)−i_(bc)=i_(b), i_(bc)−i_(ca)=i_(c) are used. From the above Equations, the differences of induced voltages of the primary windings can be obtained as follows: $\begin{matrix} {{e_{CA} - e_{AB}} = {v_{CA} - v_{AB} + {r_{1}i_{A}} + {L_{l1}\frac{i_{A}}{t}}}} & \left\lbrack {{Equation}\quad 53} \right\rbrack \\ {{e_{AB} - e_{BC}} = {v_{AB} - v_{BC} + {r_{1}i_{B}} + {L_{l1}\frac{i_{B}}{t}}}} & \left\lbrack {{Equation}\quad 54} \right\rbrack \\ {{e_{BC} - e_{CA}} = {v_{BC} - v_{CA} + {r_{1}i_{C}} + {L_{l1}\frac{i_{C}}{t}}}} & \left\lbrack {{Equation}\quad 55} \right\rbrack \end{matrix}$

Also, the differences of induced voltages of the secondary windings can be given by the following Equations. $\begin{matrix} {{e_{ca} - e_{ab}} = {v_{ca} - v_{ab} - {r_{2}i_{a}} - {L_{l2}\frac{i_{a}}{t}}}} & \left\lbrack {{Equation}\quad 56} \right\rbrack \\ {{e_{ab} - e_{bc}} = {v_{ab} - v_{bc} - {r_{2}i_{b}} - {L_{l2}\frac{i_{b}}{t}}}} & \left\lbrack {{Equation}\quad 57} \right\rbrack \\ {{e_{bc} - e_{ca}} = {v_{bc} - v_{ca} - {r_{2}i_{c}} - {L_{l2}\frac{i_{c}}{t}}}} & \left\lbrack {{Equation}\quad 58} \right\rbrack \end{matrix}$

Like Y-Δ and Δ-Y transformers, Equation 59 is a necessary and sufficient condition of Equation 60. From the above fact, the three Detectors of Equation 61, Equation 62, and Equation 63 can be used as the predetermined decision parameter. $\begin{matrix} {{\frac{e_{CA}}{e_{ca}} = \frac{N_{1}}{N_{2}}},{\frac{e_{AB}}{e_{ab}} = \frac{N_{1}}{N_{2}}},{\frac{e_{BC}}{e_{bc}} = \frac{N_{1}}{N_{2}}}} & \left\lbrack {{Equation}\quad 59} \right\rbrack \\ {{\frac{e_{CA} - e_{AB}}{e_{ca} - e_{ab}} = \frac{N_{1}}{N_{2}}},{\frac{e_{AB} - e_{BC}}{e_{ab} - e_{bc}} = \frac{N_{1}}{N_{2}}},{\frac{e_{BC} - e_{CA}}{e_{bc} - e_{ca}} = \frac{N_{1}}{N_{2}}}} & \left\lbrack {{Equation}\quad 60} \right\rbrack \\ {{{Detector}\quad 1} = {\frac{e_{CA} - e_{AB} - {\frac{N_{1}}{N_{2}}\left( {e_{ca} - e_{ab}} \right)}}{V_{CArms}} \times 100\quad (\%)}} & \left\lbrack {{Equation}\quad 61} \right\rbrack \\ {{{Detector}\quad 2} = {\frac{e_{AB} - e_{BC} - {\frac{N_{1}}{N_{2}}\left( {e_{ab} - e_{bc}} \right)}}{V_{ABrms}} \times 100\quad (\%)}} & \left\lbrack {{Equation}\quad 62} \right\rbrack \\ {{{Detector}\quad 3} = {\frac{e_{BC} - e_{CA} - {\frac{N_{1}}{N_{2}}\left( {e_{bc} - e_{ca}} \right)}}{V_{BCrms}} \times 100\quad (\%)}} & \left\lbrack {{Equation}\quad 63} \right\rbrack \end{matrix}$

The second preferred embodiment of the invention is about a protective relaying algorithm for power transformers using the ratio of flux linkage increments of primary and secondary windings. If there is no internal fault of the two-winding single-phase transformer, the primary and secondary voltages at any instant can be given by the following Equation 64 and Equation 65. $\begin{matrix} {v_{1} = {{R_{1}i_{1}} + {L_{l1}\frac{i_{1}}{t}} + \frac{\lambda_{1}}{t}}} & \left\lbrack {{Equation}\quad 64} \right\rbrack \\ {v_{2} = {{{- R_{2}}i_{2}} - {L_{l2}\frac{i_{2}}{t}} + \frac{\lambda_{2}}{t}}} & \left\lbrack {{Equation}\quad 65} \right\rbrack \end{matrix}$

where, R₁, R₂ are the primary and secondary winding resistances, respectively, and L_(l1), L_(l2) are the primary and secondary leakage inductances, respectively, and λ₁, λ₂ are primary, secondary flux linkages, respectively.

Rearranging Equation 64 and Equation 65 yields Equation 66 and Equation 67. $\begin{matrix} {\frac{\lambda_{1}}{t} = {v_{1} - {R_{1}i_{1}} - {L_{l1}\frac{i_{1}}{t}}}} & \left\lbrack {{Equation}\quad 66} \right\rbrack \\ {\frac{\lambda_{2}}{t} = {v_{2} + {R_{2}i_{2}} + {L_{l2}\frac{i_{2}}{t}}}} & \left\lbrack {{Equation}\quad 67} \right\rbrack \end{matrix}$

Integrating both sides of Equation 66 and Equation 67 from t_(n) to t_(n−1) gives the following Equation 68 and Equation 69. $\begin{matrix} {{\Delta\lambda}_{1} = {{{\lambda_{1}\left( t_{n} \right)} - {\lambda_{1}\left( t_{n - 1} \right)}} = {{\int_{n - 1}^{t_{n}}{v_{1}\quad {t}}} - {R_{1}{\int_{n - 1}^{t_{n}}{i_{1}\quad {t}}}} - {L_{l1}\left( {{i_{1}\left( t_{n} \right)} - {i_{1}\left( t_{n - 1} \right)}} \right)}}}} & \left\lbrack {{Equation}\quad 68} \right\rbrack \\ {{\Delta\lambda}_{2} = {{{\lambda_{2}\left( t_{n} \right)} - {\lambda_{2}\left( t_{n - 1} \right)}} = {{\int_{n - 1}^{t_{n}}{v_{2}\quad {t}}} + {R_{2}{\int_{n - 1}^{t_{n}}{i_{2}\quad {t}}}} - {L_{l2}\left( {{i_{2}\left( t_{n} \right)} - {i_{2}\left( t_{n - 1} \right)}} \right)}}}} & \left\lbrack {{Equation}\quad 69} \right\rbrack \end{matrix}$

Thus, the increments of the flux linkage of the primary and secondary windings Δλ₁, Δλ₂ can be estimated with voltages, currents, winding resistances and leakage inductances. The ratio of Equation 68 to Equation 69 is defined by the ratio of flux linkage increments (RFLI) represented by Equation 70. $\begin{matrix} {{RFLI} = \frac{{\Delta\lambda}_{1}}{{\Delta\lambda}_{2}}} & \left\lbrack {{Equation}\quad 70} \right\rbrack \end{matrix}$

As Δλ₁, Δλ₂ can be estimated using Equation 68 and Equation 69 the RFLI can be calculated at any instant using Equation 5. In the steady state, if there is no internal winding fault, the RFLI at any instant is equal to the turn ratio, i.e. N₁/N₂ except for Δλ₁=0 or Δλ₂=0.

In case of magnetic inrush, as the core repeats saturation and unsaturation, a magnetizing current also repeats large and small. This phenomenon lasts for some time. However, in this case, since it is not an internal winding fault, the RFLI is also equal to the turn ratio even when Δλ₁, and Δλ₂ are not sinusoidal and distorted.

On the other hand, the RFLI is not the same as the turn ratio in case of internal winding faults. Thus, from this point, the second preferred embodiment of this invention detects an internal fault of a transformer based on the following fact; it is not an internal fault if the RFLI is the same as the turn ratio, while otherwise an internal fault.

However, Δλ₁ and Δλ₂ are instantaneous values and the RFLI can grow up to a large value when Δλ₁ has values near zero. In the second preferred embodiment of this invention, a Detector of Equation 71 is used to detect a fault instead of Equation 70. That is, if Equation 71 is less than the threshold, it is normal; otherwise, it is an internal winding fault. $\begin{matrix} {{Detector} = {\frac{{\Delta\lambda}_{1} - {\frac{N_{1}}{N_{2}}{\Delta\lambda}_{2}}}{\frac{V_{1{rm}\quad s}}{4.44f}\Delta \quad T} \times 100\quad (\%)}} & \left\lbrack {{Equation}\quad 71} \right\rbrack \end{matrix}$

where, V_(1rms) is the primary rated voltage and f is nominal frequency and ΔT=t_(n)−t_(n−1).

Now, the relaying method to a two-winding three-phase Y—Y transformer is applied. In FIG. 2, voltages of the primary windings can be represented as Equation 72, Equation 73, and Equation 74. $\begin{matrix} {v_{A} = {{R_{A}i_{A}} + {L_{lA}\frac{i_{A}}{t}} + \frac{\lambda_{A}}{t}}} & \left\lbrack {{Equation}\quad 72} \right\rbrack \\ {v_{B} = {{R_{B}i_{B}} + {L_{lB}\frac{i_{B}}{t}} + \frac{\lambda_{B}}{t}}} & \left\lbrack {{Equation}\quad 73} \right\rbrack \\ {v_{C} = {{R_{C}i_{C}} + {L_{lC}\frac{i_{C}}{t}} + \frac{\lambda_{C}}{t}}} & \left\lbrack {{Equation}\quad 74} \right\rbrack \end{matrix}$

where, R_(A), R_(B), R_(C) mean primary winding resistances and L_(lA), L_(lB), L_(lC) mean primary leakage inductances, and λ_(A), λ_(B), λ_(C) mean primary flux linkages.

And, the voltages of the secondary windings are given by Equation 75, Equation 76, and Equation 77. $\begin{matrix} {v_{a} = {{{- R_{a}}i_{a}} - {L_{la}\frac{i_{a}}{t}} + \frac{\lambda_{a}}{t}}} & \left\lbrack {{Equation}\quad 75} \right\rbrack \\ {v_{b} = {{{- R_{b}}i_{b}} - {L_{l\quad b}\frac{i_{b}}{t}} + \frac{\lambda_{b}}{t}}} & \left\lbrack {{Equation}\quad 76} \right\rbrack \\ {v_{c} = {{{- R_{c}}i_{c}} - {L_{lc}\frac{i_{c}}{t}} + \frac{\lambda_{c}}{t}}} & \left\lbrack {{Equation}\quad 77} \right\rbrack \end{matrix}$

where, R_(a), R_(b), R_(c) are secondary winding resistances, and L_(la), L_(lb), L_(lc) are secondary leakage inductances, and λ_(a), λ_(b), λ_(c) are secondary flux linkages.

In case of a three-phase Y—Y transformer, in exactly the same manner of a single-phase transformer, three Detectors are defined by Equation 78, Equation 79, and Equation 80. $\begin{matrix} {{{Detector}\quad A} = {\frac{{\Delta\lambda}_{A} - {\frac{N_{1}}{N_{2}}{\Delta\lambda}_{A}}}{\frac{V_{Arms}}{4.44f}\Delta \quad T} \times 100\quad (\%)}} & \left\lbrack {{Equation}\quad 78} \right\rbrack \\ {{{Detector}\quad B} = {\frac{{\Delta\lambda}_{B} - {\frac{N_{1}}{N_{2}}{\Delta\lambda}_{b}}}{\frac{V_{Brms}}{4.44f}\Delta \quad T} \times 100\quad (\%)}} & \left\lbrack {{Equation}\quad 79} \right\rbrack \\ {{{Detector}\quad C} = {\frac{{\Delta\lambda}_{C} - {\frac{N_{1}}{N_{2}}{\Delta\lambda}_{c}}}{\frac{V_{Crms}}{4.44f}\Delta \quad T} \times 100\quad (\%)}} & \left\lbrack {{Equation}\quad 80} \right\rbrack \end{matrix}$

If values of Detector A, Detector B, and Detector C are larger than the threshold, it is determined that an internal fault occurred on phase A, phase B and phase C, respectively.

The next topic is about a two-winding three-phase Y-Δ transformer.

Voltages of the wye-connected windings (as shown in FIG. 3) are represented in Equation 72, Equation 73 and Equation 74. Equation 81, Equation 82 and Equation 83 represent voltages of delta-connected windings. $\begin{matrix} {v_{ab} = {{R_{ab}i_{ab}} + {L_{lab}\frac{i_{ab}}{t}} + \frac{\lambda_{ab}}{t}}} & \left\lbrack {{Equation}\quad 81} \right\rbrack \\ {v_{bc} = {{R_{bc}i_{bc}} + {L_{lbc}\frac{i_{bc}}{t}} + \frac{\lambda_{bc}}{t}}} & \left\lbrack {{Equation}\quad 82} \right\rbrack \\ {v_{ca} = {{R_{ca}i_{ca}} + {L_{lca}\frac{i_{ca}}{t}} + \frac{\lambda_{ca}}{t}}} & \left\lbrack {{Equation}\quad 83} \right\rbrack \end{matrix}$

where, R_(ab), R_(bc), R_(ca) are secondary winding resistances, and L_(lab), L_(lbc), L_(lca) are secondary leakage inductances, and λ_(ab), λ_(bc), λ_(ca) are secondary flux linkages.

As i_(ab), i_(bc), i_(ca) are not available in case of Y-Δ connection, λ_(ab), λ_(bc), λ_(ca) cannot be calculated directly from Equation 81, Equation 82 and Equation 83. However, i_(a), i_(b), i_(c) are available. Thus, in this invention, in order to use i_(a), i_(b), i_(c) the following relations i_(ca)−i_(ab)=i_(a), i_(ab)−i_(bc)=i_(b), i_(bc)−i_(ca)=i_(c) are employed. Moreover, if there is no internal winding fault, R_(ab)≈R_(bc)≈R_(ca)=R and L_(lab)≈L_(lbc)≈L_(lca)=L_(l). Therefore, manipulating Equation 81, Equation 82 and Equation 83 gives the following Equation 84, Equation 85 and Equation 86. $\begin{matrix} {{\frac{\lambda_{ca}}{t} - \frac{\lambda_{ab}}{t}} = {v_{ca} - v_{ab} - {Ri}_{a} - {L_{l}\frac{i_{a}}{t}}}} & \left\lbrack {{Equation}\quad 84} \right\rbrack \\ {{\frac{\lambda_{ab}}{t} - \frac{\lambda_{bc}}{t}} = {v_{ab} - v_{bc} - {Ri}_{b} - {L_{l}\frac{i_{b}}{t}}}} & \left\lbrack {{Equation}\quad 85} \right\rbrack \\ {{\frac{\lambda_{bc}}{t} - \frac{\lambda_{ca}}{t}} = {v_{bc} - v_{ca} - {Ri}_{c} - {L_{l}\frac{i_{c}}{t}}}} & \left\lbrack {{Equation}\quad 86} \right\rbrack \end{matrix}$

Therefore, left-hand sides of Equation 84, Equation 85 and Equation 86 can be calculated. The following Equation 87, Equation 88, and Equation 89 for the primary winding corresponding to Equation 84, Equation 85 and Equation 86 can be derived by manipulating Equation 72, Equation 73 and Equation 74. $\begin{matrix} {{\frac{\lambda_{C}}{t} - \frac{\lambda_{A}}{t}} = {v_{C} - v_{A} - \left( {{R_{C}i_{C}} - {R_{A}i_{A}}} \right) - \left( {{L_{lC}\frac{i_{C}}{t}} - {L_{lA}\frac{i_{A}}{t}}} \right)}} & \left\lbrack {{Equation}\quad 87} \right\rbrack \\ {{\frac{\lambda_{A}}{t} - \frac{\lambda_{B}}{t}} = {v_{A} - v_{B} - \left( {{R_{A}i_{A}} - {R_{B}i_{B}}} \right) - \left( {{L_{lA}\frac{i_{A}}{t}} - {L_{lB}\frac{i_{B}}{t}}} \right)}} & \left\lbrack {{Equation}\quad 88} \right\rbrack \\ {{\frac{\lambda_{B}}{t} - \frac{\lambda_{C}}{t}} = {v_{B} - v_{C} - \left( {{R_{B}i_{B}} - {R_{C}i_{C}}} \right) - \left( {{L_{lB}\frac{i_{B}}{t}} - {L_{lC}\frac{i_{C}}{t}}} \right)}} & \left\lbrack {{Equation}\quad 89} \right\rbrack \end{matrix}$

If there is no internal fault in case of a Y-Δ transformer, the following relationships in Equation 90 are valid. $\begin{matrix} {{\frac{{\Delta\lambda}_{A}}{{\Delta\lambda}_{ab}} = \frac{N_{1}}{N_{2}}},{\frac{{\Delta\lambda}_{B}}{{\Delta\lambda}_{bc}} = \frac{N_{1}}{N_{2}}},{\frac{{\Delta\lambda}_{C}}{{\Delta\lambda}_{ca}} = \frac{N_{1}}{N_{2}}}} & \left\lbrack {{Equation}\quad 90} \right\rbrack \end{matrix}$

While Δλ_(A), Δλ_(B), Δλ_(C) can be calculated, Δλ_(ab), Δλ_(bc), Δλ_(ca) cannot. Thus, Equation 90 cannot be applied directly to Detectors (different from the case of a Y—Y transformer). Equation 91 is used for fault detection in case of a Y-Δ transformer. $\begin{matrix} {{\frac{\Delta \left( {\lambda_{C} - \lambda_{A}} \right)}{\Delta \left( {e_{ca} - e_{ab}} \right)} = \frac{N_{1}}{N_{2}}},{\frac{\Delta \left( {\lambda_{A} - \lambda_{B}} \right)}{\Delta \left( {e_{ab} - e_{bc}} \right)} = \frac{N_{1}}{N_{2}}},{\frac{\Delta \left( {\lambda_{B} - \lambda_{C}} \right)}{\Delta \left( {e_{bc} - e_{ca}} \right)} = \frac{N_{1}}{N_{2}}}} & \left\lbrack {{Equation}\quad 91} \right\rbrack \end{matrix}$

Equation 90 is a necessary and sufficient condition of Equation 91. Sufficient condition is proved trivially in substituting Equation 90 into Equation 91. Thus, only necessary condition will be proved. Let three RFLI as Δλ_(A)/Δλ_(ab)=α, Δλ_(B)/Δλ_(bc)=β, and Δλ_(C)/Δλ_(ca)=γ. Substituting the three RFLIs into Equation 91 yields the following three equations Δ(γλ_(ca)−αλ_(ab))=(N₁/N₂)Δ(λ_(ca)−λ_(ab)), Δ(αλ_(ab)−βλ_(bc))=(N₁/N₂)Δ(λ_(ab)−λ_(bc)), and Δ(βλ_(bc)−γλ_(ca))=(N₁/N₂)Δ(λ_(bc)−λ_(ca)). In order for the former three Equations to be valid for all Δλ_(ab), Δλ_(bc), Δλ_(ca) the condition of α=β=γ=N₁/N₂ should be satisfied. Thus, Δλ_(A)/Δλ_(ab)=Δλ_(B)/Δλ_(bc)=Δλ_(C)/Δλ_(ca)=N₁/N₂. Therefore, Equation 90 is a necessary and sufficient condition of Equation 91.

Thus, in the second preferred embodiment of the invention, the Detectors for a Y-Δ transformer are given by Equation 92, Equation 93, and Equation 94. $\begin{matrix} {{{Detector}\quad 1} = {\frac{{\Delta \left( {\lambda_{C} - \lambda_{A}} \right)} - {\frac{N_{1}}{N_{2}}{\Delta \left( {\lambda_{ca} - \lambda_{ab}} \right)}}}{\frac{V_{CArms}}{4.44f}\Delta \quad T} \times 100\quad (\%)}} & \left\lbrack {{Equation}\quad 92} \right\rbrack \\ {{{Detector}\quad 2} = {\frac{{\Delta \left( {\lambda_{A} - \lambda_{B}} \right)} - {\frac{N_{1}}{N_{2}}{\Delta \left( {\lambda_{ab} - \lambda_{bc}} \right)}}}{\frac{V_{ABrms}}{4.44f}\Delta \quad T} \times 100\quad (\%)}} & \left\lbrack {{Equation}\quad 93} \right\rbrack \end{matrix}$

$\begin{matrix} {{{Detector}\quad 3} = {\frac{{\Delta \left( {\lambda_{B} - \lambda_{C}} \right)} - {\frac{N_{1}}{N_{2}}{\Delta \left( {\lambda_{bc} - \lambda_{ca}} \right)}}}{\frac{V_{BCrms}}{4.44f}\Delta \quad T} \times 100\quad (\%)}} & \left\lbrack {{Equation}\quad 94} \right\rbrack \end{matrix}$

The rule of detecting internal faults and the faulted phase from the above three Detectors in case of a Y-Δ transformer is the same as Table 1 and its detailed explanation is omitted.

Case studies of the second preferred embodiment of a protective relaying method conducted on the model system of FIG. 4 is discussed.

The first case is magnetic inrush.

We will explain the two magnetic inrush cases of with 0(%) and +80(%) remanent fluxes of the saturation point in case of 0 (deg) energization angle and no load. FIG. 5 shows three-phase currents in case of magnetic inrush with a 0(%) remanent flux, 0 (deg) energization angle, and no load.

FIG. 11 shows calculated flux linkage difference's increments of Detector 1 i.e. Δ(λ_(C)−λ_(A)) and (N₁/N₂)Δ(λ_(ca)−λ_(ab)), increments of Detector 2 i.e. Δ(λ_(A)−λ_(B)) and (N₁/N₂)Δ(λ_(ab)−λ_(bc)), and increments of Detector 3 i.e. Δ(λ_(B)−λ_(C)) and (N₁/N₂)Δ(λ_(bc)−λ_(ca)). The result indicates that the increments of the calculated flux linkage difference are nearly the same even if inrush currents contain harmonic components.

FIG. 12 shows three Detectors and a trip signal in case of no remanent flux according to the second preferred embodiment of this invention. The trip signal in FIG. 12 means the signal to initiate the circuit breaker after detection of a fault.

The fault detection rule is the same as that of the first case. As the three Detectors are less than 5(%), the trip signal is blocked even the A-phase inrush current increases up to nearly 900 (A) after the energization.

The errors shown in Detectors of FIG. 12 are small, similar to those of the first embodiment.

FIG. 13 and FIG. 14 show currents of each phase, the three Detectors, and a trip signal in case of +80% remanent flux of the saturation point in the transformer. Like the first preferred embodiment, all the three Detectors are less than 5(%) even in case of very large inrush currents and the trip signal is blocked.

Internal winding fault is now explained.

FIG. 15 shows the three Detectors and a trip signal in case of the B-phase internal turn-to-earth fault at the 40% point from the neutral according to the second preferred embodiment of this invention. Like the first embodiment, the proposed method detects the internal winding fault 4.167 (ms) after a fault. Finally, the second preferred embodiment proposes a transformer protective relaying algorithm using the ratio of flux linkage increments of the primary and secondary windings. The proposed algorithm calculates the flux linkage increments from the currents and voltages and then clearly discriminates between internal winding faults and magnetic inrush with the ratio. Various test results clearly demonstrate that the algorithm does not malfunction even in case of inrush with the maximum remanent flux. Moreover, it detects internal winding faults in only ¼ cycle after a fault.

Although the invention is described herein with reference to the preferred embodiment, it should be clearly understood that many variations and/or modifications of the basic inventive concepts herein taught which may be apparent to those skilled in the present art will remain within the spirit and scope of the present invention, as defined in the appended claims.

The cases of three-phase Δ-Y transformer and Δ-Δ transformer are briefly described below as examples of the above mentioned modifications.

FIG. 17 shows a two-winding three-phase Δ-Y transformer. Voltages of the primary windings can be represented by the following Equations. $\begin{matrix} {v_{AB} = {{r_{AB}i_{AB}} + {L_{lAB}\frac{i_{AB}}{t}} + \frac{\lambda_{AB}}{t}}} & \left\lbrack {{Equation}\quad 95} \right\rbrack \\ {v_{BC} = {{r_{BC}i_{BC}} + {L_{lBC}\frac{i_{BC}}{t}} + \frac{\lambda_{BC}}{t}}} & \left\lbrack {{Equation}\quad 96} \right\rbrack \\ {v_{CA} = {{r_{CA}i_{CA}} + {L_{lCA}\frac{i_{CA}}{t}} + \frac{\lambda_{cA}}{t}}} & \left\lbrack {{Equation}\quad 97} \right\rbrack \end{matrix}$

Also, voltages of the secondary windings can be represented by the following Equations. $\begin{matrix} {v_{a} = {{{- R_{a}}i_{a}} - {L_{la}\frac{i_{a}}{t}} + \frac{\lambda_{a}}{t}}} & \left\lbrack {{Equation}\quad 98} \right\rbrack \\ {v_{b} = {{{- R_{b}}i_{b}} - {L_{l\quad b}\frac{i_{b}}{t}} + \frac{\lambda_{b}}{t}}} & \left\lbrack {{Equation}\quad 99} \right\rbrack \\ {v_{c} = {{{- R_{c}}i_{c}} - {L_{lc}\frac{i_{c}}{t}} + \frac{\lambda_{c}}{t}}} & \left\lbrack {{Equation}\quad 100} \right\rbrack \end{matrix}$

In a similar manner to a Y-Δ transformer, the relationships of r_(AB)≈r_(BC)≈r_(CA)≈r, L_(lAB)≈L_(lBC)≈L_(lCA)≈L_(l), i_(CA)−i_(AB)=−i_(A), i_(AB)−i_(BC)=−i_(B), i_(BC)−i_(CA)=−i_(C) are used. From the above Equations, flux tO linkage difference's increments the primary windings can be obtained as follows: $\begin{matrix} {{\frac{\lambda_{CA}}{t} - \frac{\lambda_{AB}}{t}} = {v_{CA} - v_{AB} + {ri}_{A} + {L_{l}\frac{i_{A}}{t}}}} & \left\lbrack {{Equation}\quad 101} \right\rbrack \\ {{\frac{\lambda_{AB}}{t} - \frac{\lambda_{BC}}{t}} = {v_{AB} - v_{BC} + {ri}_{B} + {L_{l}\frac{i_{B}}{t}}}} & \left\lbrack {{Equation}\quad 102} \right\rbrack \\ {{\frac{\lambda_{BC}}{t} - \frac{\lambda_{CA}}{t}} = {v_{BC} - v_{CA} + {ri}_{C} + {L_{l}\frac{i_{C}}{t}}}} & \left\lbrack {{Equation}\quad 103} \right\rbrack \end{matrix}$

Also, the differences of induced voltages of the secondary windings can be given by the following Equations. $\begin{matrix} {{\frac{\lambda_{c}}{t} - \frac{\lambda_{a}}{t}} = {v_{c} - v_{a} + \left( {{r_{c}i_{c}} - {r_{a}i_{a}}} \right) + \left( {{L_{l\quad c}\frac{i_{c}}{t}} - {L_{l\quad a}\frac{i_{a}}{t}}} \right)}} & \left\lbrack {{Equation}\quad 104} \right\rbrack \\ {{\frac{\lambda_{a}}{t} - \frac{\lambda_{b}}{t}} = {v_{a} - v_{b} + \left( {{r_{a}i_{a}} - {r_{b}i_{b}}} \right) + \left( {{L_{l\quad a}\frac{i_{a}}{t}} - {L_{l\quad b}\frac{i_{b}}{t}}} \right)}} & \left\lbrack {E\quad q\quad u\quad a\quad t\quad i\quad o\quad n\quad 105} \right\rbrack \\ {{\frac{\lambda_{b}}{t} - \frac{\lambda_{c}}{t}} = {v_{b} - v_{c} + \left( {{r_{b}i_{b}} - {r_{c}i_{c}}} \right) + \left( {{L_{l\quad b}\frac{i_{b}}{t}} - {L_{l\quad c}\frac{i_{c}}{t}}} \right)}} & \left\lbrack {E\quad q\quad u\quad a\quad t\quad i\quad o\quad n\quad 106} \right\rbrack \end{matrix}$

Like a Y-Δ transformer, Equation 107 is a necessary and sufficient condition of Equation 108. From the above fact, the three Detectors of Equation 109, Equation 110, and Equation 111 can be used as the predetermined decision parameter. $\begin{matrix} {{\frac{{\Delta\lambda}_{C\quad A}}{{\Delta\lambda}_{c}} = \frac{N_{1}}{N_{2}}},{\frac{{\Delta\lambda}_{A\quad B}}{{\Delta\lambda}_{a}} = \frac{N_{1}}{N_{2}}},{\frac{{\Delta\lambda}_{B\quad C}}{{\Delta\lambda}_{b}} = \frac{N_{1}}{N_{2}}}} & \left\lbrack {{Equation}\quad 107} \right\rbrack \\ {{\frac{{\Delta\lambda}_{C\quad A} - {\Delta\lambda}_{A\quad B}}{{\Delta\lambda}_{c} - {\Delta\lambda}_{a}} = \frac{N_{1}}{N_{2}}},{\frac{{\Delta\lambda}_{A\quad B} - {\Delta\lambda}_{B\quad C}}{{\Delta\lambda}_{a} - {\Delta\lambda}_{b}} = \frac{N_{1}}{N_{2}}},{\frac{{\Delta\lambda}_{B\quad C} - {\Delta\lambda}_{C\quad A}}{{\Delta\lambda}_{b} - {\Delta\lambda}_{c}} = \frac{N_{1}}{N_{2}}}} & \left\lbrack {E\quad q\quad u\quad a\quad t\quad i\quad o\quad n\quad 108} \right\rbrack \\ {{D\quad e\quad t\quad e\quad c\quad t\quad o\quad r\quad 1} = {\frac{{\Delta \left( {\lambda_{C\quad A} - \lambda_{A\quad B}} \right)} - {\frac{N_{1}}{N_{2}}{\Delta \left( {\lambda_{c} - \lambda_{a}} \right)}}}{\frac{V_{C\quad A\quad r\quad m\quad s}}{4.44f}\Delta \quad T} \times 100(\%)}} & \left\lbrack {E\quad q\quad u\quad a\quad t\quad i\quad o\quad n\quad 109} \right\rbrack \\ {{D\quad e\quad t\quad e\quad c\quad t\quad o\quad r\quad 2} = {\frac{{\Delta \left( {\lambda_{A\quad B} - \lambda_{B\quad C}} \right)} - {\frac{N_{1}}{N_{2}}{\Delta \left( {\lambda_{a} - \lambda_{b}} \right)}}}{\frac{V_{A\quad B\quad r\quad m\quad s}}{4.44f}\Delta \quad T} \times 100(\%)}} & \left\lbrack {E\quad q\quad u\quad a\quad t\quad i\quad o\quad n\quad 110} \right\rbrack \\ {{D\quad e\quad t\quad e\quad c\quad t\quad o\quad r\quad 3} = {\frac{{\Delta \left( {\lambda_{B\quad C} - \lambda_{C\quad A}} \right)} - {\frac{N_{1}}{N_{2}}{\Delta \left( {\lambda_{b} - \lambda_{c}} \right)}}}{\frac{V_{B\quad C\quad r\quad m\quad s}}{4.44f}\Delta \quad T} \times 100(\%)}} & \left\lbrack {E\quad q\quad u\quad a\quad t\quad i\quad o\quad n\quad 111} \right\rbrack \end{matrix}$

FIG. 16 shows a two-winding three-phase Δ-Δ transformer. Voltages of the primary windings can be represented by the following Equations. $\begin{matrix} {v_{A\quad B} = {{r_{A\quad B}i_{A\quad B}} + {L_{l\quad A\quad B}\frac{i_{A\quad B}}{t}} + \frac{\lambda_{A\quad B}}{t}}} & \left\lbrack {{Equation}\quad 112} \right\rbrack \\ {v_{B\quad C} = {{r_{B\quad C}i_{B\quad C}} + {L_{l\quad B\quad C}\frac{i_{B\quad C}}{t}} + \frac{\lambda_{B\quad C}}{t}}} & \left\lbrack {E\quad q\quad u\quad a\quad t\quad i\quad o\quad n\quad 113} \right\rbrack \\ {v_{C\quad A} = {{r_{C\quad A}i_{C\quad A}} + {L_{l\quad C\quad A}\frac{i_{C\quad A}}{t}} + \frac{\lambda_{C\quad A}}{t}}} & \left\lbrack {E\quad q\quad u\quad a\quad t\quad i\quad o\quad n\quad 114} \right\} \end{matrix}$

Also, voltages of the secondary windings can be represented by the following Equations. $\begin{matrix} {v_{a\quad b} = {{r_{a\quad b}i_{a\quad b}} + {L_{l\quad a\quad b}\frac{i_{a\quad b}}{t}} + \frac{\lambda_{a\quad b}}{t}}} & \left\lbrack {{Equation}\quad 115} \right\rbrack \\ {v_{b\quad c} = {{r_{b\quad c}i_{b\quad c}} + {L_{l\quad b\quad c}\frac{i_{b\quad c}}{t}} + \frac{\lambda_{b\quad c}}{t}}} & \left\lbrack {E\quad q\quad u\quad a\quad t\quad i\quad o\quad n\quad 116} \right\rbrack \\ {v_{c\quad a} = {{r_{c\quad a}i_{c\quad a}} + {L_{l\quad c\quad a}\frac{i_{c\quad a}}{t}} + \frac{\lambda_{c\quad a}}{t}}} & \left\lbrack {E\quad q\quad u\quad a\quad t\quad i\quad o\quad n\quad 117} \right\} \end{matrix}$

In a similar manner to a Δ-Y transformer, the relationships of r_(AB)≈r_(BC)≈r_(CA)≈r₁, L_(lAB)≈L_(lBC)≈L_(lCA)≈L_(l1), i_(CA)−i_(AB)=−i_(A), i_(AB)−i_(BC)=−i_(B), i_(BC)−i_(CA)=−i_(C), r_(ab)≈r_(bc)≈r_(ca)≈r₂, L_(lab)≈L_(lbc)≈L_(lca)≈L_(l2), i_(ca)−i_(ab)=i_(a), i_(ab)−i_(bc)=i_(b), i_(bc)−i_(ca)=i_(c) are used. From the above Equations, the flux linkage difference's increments of the primary windings can be obtained as follows: $\begin{matrix} {{\frac{\lambda_{C\quad A}}{t} - \frac{\lambda_{A\quad B}}{t}} = {v_{C\quad A} - v_{A\quad B} + {r_{1}i_{A}} + {L_{l1}\frac{i_{A}}{t}}}} & \left\lbrack {{Equation}\quad 118} \right\rbrack \\ {{\frac{\lambda_{A\quad B}}{t} - \frac{\lambda_{B\quad C}}{t}} = {v_{A\quad B} - v_{B\quad C} + {r_{1}i_{B}} + {L_{l1}\frac{i_{B}}{t}}}} & \left\lbrack {E\quad q\quad u\quad a\quad t\quad i\quad o\quad n\quad 119} \right\rbrack \\ {{\frac{\lambda_{B\quad C}}{t} - \frac{\lambda_{C\quad A}}{t}} = {v_{B\quad C} - v_{C\quad A} + {r_{1}i_{C}} + {L_{l1}\frac{i_{C}}{t}}}} & \left\lbrack {E\quad q\quad u\quad a\quad t\quad i\quad o\quad n\quad 120} \right\rbrack \end{matrix}$

Also, the differences of induced voltages of the secondary windings can be given by the following Equations. $\begin{matrix} {{\frac{\lambda_{c\quad a}}{t} - \frac{\lambda_{a\quad b}}{t}} = {v_{c\quad a} - v_{a\quad b} - {r_{2}i_{a}} - {L_{l2}\frac{i_{a}}{t}}}} & \left\lbrack {{Equation}\quad 121} \right\rbrack \\ {{\frac{\lambda_{a\quad b}}{t} - \frac{\lambda_{b\quad c}}{t}} = {v_{a\quad b} - v_{b\quad c} - {r_{2}i_{b}} - {L_{l2}\frac{i_{b}}{t}}}} & \left\lbrack {E\quad q\quad u\quad a\quad t\quad i\quad o\quad n\quad 122} \right\rbrack \\ {{\frac{\lambda_{b\quad c}}{t} - \frac{\lambda_{c\quad a}}{t}} = {v_{b\quad c} - v_{c\quad a} - {r_{2}i_{c}} - {L_{l2}\frac{i_{c}}{t}}}} & \left\lbrack {E\quad q\quad u\quad a\quad t\quad i\quad o\quad n\quad 123} \right\rbrack \end{matrix}$

Like Y-Δ and Δ-Y transformers, Equation 124 is a necessary and sufficient condition of Equation 125. From the above fact, the three Detectors of Equation 126, Equation 127, and Equation 128 can be used as the predetermined decision parameter. $\begin{matrix} {{\frac{{\Delta\lambda}_{C\quad A}}{{\Delta\lambda}_{c\quad a}} = \frac{N_{1}}{N_{2}}},{\frac{{\Delta\lambda}_{A\quad B}}{{\Delta\lambda}_{a\quad b}} = \frac{N_{1}}{N_{2}}},{\frac{{\Delta\lambda}_{B\quad C}}{{\Delta\lambda}_{b\quad c}} = \frac{N_{1}}{N_{2}}}} & \left\lbrack {{Equation}\quad 124} \right\rbrack \\ {{\frac{{\Delta\lambda}_{C\quad A} - {\Delta\lambda}_{A\quad B}}{{\Delta\lambda}_{c\quad a} - {\Delta\lambda}_{a\quad b}} = \frac{N_{1}}{N_{2}}},{\frac{{\Delta\lambda}_{A\quad B} - {\Delta\lambda}_{B\quad C}}{{\Delta\lambda}_{a\quad b} - {\Delta\lambda}_{b\quad c}} = \frac{N_{1}}{N_{2}}},{\frac{{\Delta\lambda}_{B\quad C} - {\Delta\lambda}_{C\quad A}}{{\Delta\lambda}_{b\quad c} - {\Delta\lambda}_{c\quad a}} = \frac{N_{1}}{N_{2}}}} & \left\lbrack {E\quad q\quad u\quad a\quad t\quad i\quad o\quad n\quad 125} \right\rbrack \\ {{D\quad e\quad t\quad e\quad c\quad t\quad o\quad r\quad 1} = {\frac{{\Delta \left( {\lambda_{C\quad A} - \lambda_{A\quad B}} \right)} - {\frac{N_{1}}{N_{2}}{\Delta \left( {\lambda_{c\quad a} - \lambda_{a\quad b}} \right)}}}{\frac{V_{C\quad A\quad r\quad m\quad s}}{4.44f}\Delta \quad T} \times 100(\%)}} & \left\lbrack {E\quad q\quad u\quad a\quad t\quad i\quad o\quad n\quad 1\quad 26} \right\rbrack \\ {{D\quad e\quad t\quad e\quad c\quad t\quad o\quad r\quad 2} = {\frac{{\Delta \left( {\lambda_{A\quad B} - \lambda_{B\quad C}} \right)} - {\frac{N_{1}}{N_{2}}{\Delta \left( {\lambda_{a\quad b} - \lambda_{b\quad c}} \right)}}}{\frac{V_{A\quad B\quad r\quad m\quad s}}{4.44f}\Delta \quad T} \times 100(\%)}} & \left\lbrack {E\quad q\quad u\quad a\quad t\quad i\quad o\quad n\quad 127} \right\rbrack \\ {{D\quad e\quad t\quad e\quad c\quad t\quad o\quad r\quad 3} = {\frac{{\Delta \left( {\lambda_{B\quad C} - \lambda_{C\quad A}} \right)} - {\frac{N_{1}}{N_{2}}{\Delta \left( {\lambda_{b\quad c} - \lambda_{c\quad a}} \right)}}}{\frac{V_{B\quad C\quad r\quad m\quad s}}{4.44f}\Delta \quad T} \times 100(\%)}} & \left\lbrack {E\quad q\quad u\quad a\quad t\quad i\quad o\quad n\quad 128} \right\rbrack \end{matrix}$

Although the invention is described herein with reference to the preferred embodiment, it should be clearly understood that many variations and/or modifications of the basic inventive concepts herein taught which may be apparent to those skilled in the present art will remain within the spirit and scope of the present invention, as defined in the appended claims.

The invention proposes a transformer protective relaying algorithm using the ratio of induced voltages or ratio of flux linkage increments. The ratios are the same as the turn ratio in case of magnetic inrush and overexcitation while they are not equal to the turn ratio in case of internal winding faults. Thus, The proposed algorithm calculates the induced voltages and flux linkage increments of the primary and secondary windings from the currents and voltages and then discriminates between internal winding faults and magnetic inrush with the ratio.

The proposed method does not require hysteresis data and can be applied without significant modification even if system and core conditions are significantly changed. In addition to that, as the algorithm works in time domain, it can reduce the relay's operating time, especially in a higher sampling rate. 

What is claimed is:
 1. A relaying method for protecting a transformer comprising the steps of: obtaining primary and secondary and voltages and currents of said transformer; calculating induced voltages, induced voltage differences, ratio of primary and secondary induced voltages, or ratio of primary and secondary induced voltage differences from said currents and said voltages; calculating at least one predetermined decision parameter derived from at least one predetermined eguation; deciding whether an internal winding fault occurs by comparing said decision parameter to said induced voltages, said induced voltage differences, said ratio of primary and secondary induced voltages, or said ratio of primary and secondary induced voltage differences; wherein said transformer is a single-phase transformer; and wherein, in the step of calculating at least one predetermined decision parameter derived from at least one predetermined equation, said predetermined equation is $\frac{e_{1}}{e_{2}} = \frac{N_{1}}{N_{2}}$

or a modification of said equation (e₁, e₂: primary and secondary induced voltages; N₁, N₂: numbers of turns of primary and secondary windings).
 2. The method of claim 1, wherein said decision parameter is: $\frac{e_{1} - {\frac{N_{1}}{N_{2}}e_{2}}}{V_{1{rms}}} \times 100\% {\left( {V_{1{rms}}\text{:}{primary}\quad {rated}\quad {voltage}} \right).}$


3. A relaying method for protecting a transformer comprising the steps of: obtaining primary and secondary and voltages and currents of said transformer; calculating induced voltages, induced voltage differences, ratio of primary and secondary induced voltages, or ratio of primary and secondary induced voltage differences from said currents and said voltages; calculating at least one predetermined decision parameter derived from at least one predetermined equation; deciding whether an internal winding fault occurs by comparing said decision parameter to said induced voltages, said induced voltage differences, said ratio of primary and secondary induced voltages, or said ratio of primary and secondary induced voltage differences; wherein said transformer is a three-phase Y—Y transformer; and wherein, in the step of calculating at least one predetermined decision parameter derived from at least one predetermined equation, said predetermined equations are ${\frac{e_{A}}{e_{a}} = \frac{N_{1}}{N_{2}}},\quad {\frac{e_{B}}{e_{b}} = \frac{N_{1}}{N_{2}}},\quad {\frac{e_{C}}{e_{c}} = \frac{N_{1}}{N_{2}}}$

or modifications of said equations (e_(A), e_(B), e_(c): primary induced voltages; e_(a), e_(b), e_(c): secondary induced voltages; N₁, N₂: numbers of turns of primary and secondary windings).
 4. The method of claim 3, wherein said decision parameters are: ${{Detector}\quad A} = {\frac{e_{A} - {\frac{N_{1}}{N_{2}}e_{a}}}{V_{Arms}} \times 100(\%)}$ ${{Detector}\quad B} = {\frac{e_{B} - {\frac{N_{1}}{N_{2}}e_{b}}}{V_{Brms}} \times 100(\%)}$ ${{Detector}\quad C} = {\frac{e_{C} - {\frac{N_{1}}{N_{2}}e_{c}}}{V_{Crms}} \times 100(\%)}$ (V_(Arms), V_(Brms), V_(Crms):primary  rated  voltages).


5. A relaying method for protecting a transformer comprising the steps of: obtaining primary and secondary and voltages and currents of said transformer; calculating induced voltages, induced voltage differences, ratio of primary and secondary induced voltages, or ratio of primary and secondary induced voltage differences from said currents and said voltages; calculating at least one predetermined decision parameter derived from at least one predetermined equation; deciding whether an internal winding fault occurs by comparing said decision parameter to said induced voltages, said induced voltage differences, said ratio of primary and secondary induced voltages, or said ratio of primary and secondary induced voltage differences; wherein said transformer is a three-phase Y-Δ transformer; and wherein, in the step of calculating at least one predetermined decision parameter derived from at least one predetermined equation, said predetermined equations are ${\frac{e_{A}}{e_{ab}} = \frac{N_{1}}{N_{2}}},\quad {\frac{e_{B}}{e_{bc}} = \frac{N_{1}}{N_{2}}},\quad {\frac{e_{C}}{e_{ca}} = {{\frac{N_{1}}{N_{2}}\quad {or}\quad \frac{e_{C} - e_{A}}{e_{ca} - e_{ab}}} = \frac{N_{1}}{N_{2}}}},\quad {\frac{e_{A} - e_{B}}{e_{ab} - e_{bc}} = \frac{N_{1}}{N_{2}}},\quad {\frac{e_{B} - e_{C}}{e_{bc} - e_{ca}} = \frac{N_{1}}{N_{2}}}$

or modifications of said equations (e_(A), e_(B), e_(C): primary induced voltages; e_(ab), e_(bc), e_(ca): secondary induced voltages; N₁, N₂: numbers of turns of primary and secondary windings).
 6. The method of claim 5, wherein said decision parameters are: ${{Detector}\quad 1} = {\frac{e_{C} - e_{A} - {\frac{N_{1}}{N_{2}}\left( {e_{ca} - e_{ab}} \right)}}{V_{CArms}} \times 100(\%)}$ ${{Detector}\quad 2} = {\frac{e_{A} - e_{B} - {\frac{N_{1}}{N_{2}}\left( {e_{ab} - e_{bc}} \right)}}{V_{ABrms}} \times 100(\%)}$ ${{Detector}\quad 3} = {\frac{e_{B} - e_{C} - {\frac{N_{1}}{N_{2}}\left( {e_{bc} - e_{ca}} \right)}}{V_{BCrms}} \times 100(\%)}$ (V_(ABrms), V_(BCrms), V_(CArms):primary  rated  voltages).


7. A relaying method for protecting a transformer comprising the steps of: obtaining primary and secondary and voltages and currents of said transformer; calculating induced voltages, induced voltage differences, ratio of primary and secondary induced voltages, or ratio of primary and secondary induced voltage differences from said currents and said voltages; calculating at least one predetermined decision parameter derived from at least one predetermined equation; deciding whether an internal winding fault occurs by comparing said decision parameter to said induced voltages, said induced voltage differences, said ratio of primary and secondary induced voltages, or said ratio of primary and secondary induced voltage differences; wherein said transformer is a three-phase Δ-Y transformer; and wherein, in the step of calculating at least one predetermined decision parameter derived from at least one predetermined equation, said predetermined equations are ${\frac{e_{CA}}{e_{c}} = \frac{N_{1}}{N_{2}}},\quad {\frac{e_{AB}}{e_{a}} = \frac{N_{1}}{N_{2}}},\quad {\frac{e_{BC}}{e_{b}} = {{\frac{N_{1}}{N_{2}}\quad {or}\quad \frac{e_{CA} - e_{AB}}{e_{c} - e_{a}}} = \frac{N_{1}}{N_{2}}}},\quad {\frac{e_{AB} - e_{BC}}{e_{a} - e_{b}} = \frac{N_{1}}{N_{2}}},\quad {\frac{e_{BC} - e_{CA}}{e_{b} - e_{c}} = \frac{N_{1}}{N_{2}}}$

or modifications of said equations (e_(a), e_(b), e_(c): secondary induced voltages; e_(AB), e_(BC), e_(CA): primary induced voltages; N₁, N₂: numbers of turns of primary and secondary windings).
 8. The method of claim 7, wherein said decision parameters are: ${{Detector}\quad 1} = {\frac{e_{CA} - e_{AB} - {\frac{N_{1}}{N_{2}}\left( {e_{c} - e_{a}} \right)}}{V_{CArms}} \times 100(\%)}$ ${{Detector}\quad 2} = {\frac{e_{AB} - e_{BC} - {\frac{N_{1}}{N_{2}}\left( {e_{a} - e_{b}} \right)}}{V_{ABrms}} \times 100(\%)}$ ${{Detector}\quad 3} = {\frac{e_{BC} - e_{CA} - {\frac{N_{1}}{N_{2}}\left( {e_{b} - e_{c}} \right)}}{V_{BCrms}} \times 100(\%)}$ (V_(ABrms), V_(BCrms), V_(CArms):primary  rated  voltages).


9. A relaying method for protecting a transformer comprising the steps of: obtaining primary and secondary and voltages and currents of said transformer; calculating induced voltages, induced voltage differences, ratio of primary and secondary induced voltages, or ratio of primary and secondary induced voltage differences from said currents and said voltages; calculating at least one predetermined decision parameter derived from at least one predetermined equation; deciding whether an internal winding fault occurs by comparing said decision parameter to said induced voltages, said induced voltage differences, said ratio of primary and secondary induced voltages, or said ratio of primary and secondary induced voltage differences; wherein said transformer is a three-phase Δ-Δ transformer; and wherein, in the step of calculating at least one predetermined decision parameter derived from at least one predetermined equation, said predetermined equations are ${\frac{e_{CA}}{e_{ca}} = \frac{N_{1}}{N_{2}}},\quad {\frac{e_{AB}}{e_{ab}} = \frac{N_{1}}{N_{2}}},\quad {\frac{e_{BC}}{e_{bc}} = {{\frac{N_{1}}{N_{2}}\quad {or}\quad \frac{e_{CA} - e_{AB}}{e_{ca} - e_{ab}}} = \frac{N_{1}}{N_{2}}}},\quad {\frac{e_{AB} - e_{BC}}{e_{ab} - e_{bc}} = \frac{N_{1}}{N_{2}}},\quad {\frac{e_{BC} - e_{CA}}{e_{bc} - e_{ca}} = \frac{N_{1}}{N_{2}}}$

or modifications of said equations (e_(ab), e_(bc), e_(ca): secondary induced voltages; e_(AB), e_(BC), e_(CA): primary induced voltages; N₁, N₂: numbers of turns of primary and secondary windings).
 10. The method of claim 9, wherein said decision parameters are: ${{Detector}\quad 1} = {\frac{e_{CA} - e_{AB} - {\frac{N_{1}}{N_{2}}\left( {e_{ca} - e_{ab}} \right)}}{V_{CArms}} \times 100(\%)}$ ${{Detector}\quad 2} = {\frac{e_{AB} - e_{BC} - {\frac{N_{1}}{N_{2}}\left( {e_{ab} - e_{bc}} \right)}}{V_{ABrms}} \times 100(\%)}$ ${{Detector}\quad 3} = {\frac{e_{BC} - e_{CA} - {\frac{N_{1}}{N_{2}}\left( {e_{bc} - e_{ca}} \right)}}{V_{BCrms}} \times 100(\%)}$ (V_(ABrms), V_(BCrms), V_(CArms):primary  rated  voltages).


11. A relaying method for protecting a transformer comprising the steps of: obtaining primary and secondary and voltages and currents of said transformer; calculating flux linkage increments, flux linkage difference increments, the ratio of flux linkage increments of the primary and secondary windings, or the ratio of flux linkage difference increments of the primary and secondary windings from said currents and said voltages; calculating at least one predetermined decision parameter derived from at least one predetermined equation; deciding whether an internal winding fault occurs by comparing said decision parameter to said flux linkage increments, said flux linkage difference increments, said ratio of flux linkage increments of the primary and secondary windings, or said ratio of flux linkage difference increments of the primary and secondary windings; wherein said transformer is a single-phase transformer; and wherein, in the step of calculating at least one predetermined decision parameter derived from at least one predetermined equation, said predetermined equation is $\frac{{\Delta\lambda}_{1}}{{\Delta\lambda}_{2}} = \frac{N_{1}}{N_{2}}$

or a modification of said equation (Δλ₁, Δλ₂: flux linkage increments of primary and secondary windings; N₁, N₂: numbers of turns of primary and secondary windings).
 12. The method of claim 11, wherein said decision parameter is: ${Detector} = {\frac{{\Delta\lambda}_{1} - {\frac{N_{1}}{N_{2}}{\Delta\lambda}_{2}}}{\frac{V_{1{rms}}}{4.44f}\Delta \quad T} \times 100(\%){\left( {V_{1{rms}}\text{:}{primary}\quad {rated}\quad {voltage}} \right).}}$


13. A relaying method for protecting a transformer comprising the steps of: obtaining primary and secondary and voltages and currents of said transformer; calculating flux linkage increments, flux linkage difference increments, the ratio of flux linkage increments of the primary and secondary windings, or the ratio of flux linkage difference increments of the primary and secondary windings from said currents and said voltages; calculating at least one predetermined decision parameter derived from at least one predetermined equation; deciding whether an internal winding fault occurs by comparing said decision parameter to said flux linkage increments, said flux linkage difference increments, said ratio of flux linkage increments of the primary and secondary windings, or said ratio of flux linkage difference increments of the primary and secondary windings; wherein said transformer is a three-phase Y—Y transformer; and wherein, in the step of calculating at least one predetermined decision parameter derived from at least one predetermined equation, said predetermined equations are ${\frac{{\Delta\lambda}_{A}}{{\Delta\lambda}_{a}} = \frac{N_{1}}{N_{2}}},\quad {\frac{{\Delta\lambda}_{B}}{{\Delta\lambda}_{b}} = \frac{N_{1}}{N_{2}}},\quad {\frac{{\Delta\lambda}_{C}}{{\Delta\lambda}_{c}} = \frac{N_{1}}{N_{2}}}$

or modifications of said equations (Δλ_(A), Δλ_(B), Δλ_(C): primary flux linkage increments; Δλ_(a), Δλ_(b), Δλ_(c): secondary flux linkage increments; N₁, N₂: numbers of turns of primary and secondary windings).
 14. The method of claim 13, wherein said decision parameters are: ${{Detector}\quad A} = {\frac{{\Delta\lambda}_{A} - {\frac{N_{1}}{N_{2}}{\Delta\lambda}_{a}}}{\frac{V_{Arms}}{4.44f}\Delta \quad T} \times 100(\%)}$ ${{Detector}\quad B} = {\frac{{\Delta\lambda}_{B} - {\frac{N_{1}}{N_{2}}{\Delta\lambda}_{b}}}{\frac{V_{Brms}}{4.44f}\Delta \quad T} \times 100(\%)}$ ${{Detector}\quad C} = {\frac{{\Delta\lambda}_{C} - {\frac{N_{1}}{N_{2}}{\Delta\lambda}_{c}}}{\frac{V_{Crms}}{4.44f}\Delta \quad T} \times 100(\%)}$ (V_(Arms), V_(Brms), V_(Crms):primary  rated  voltages).


15. A relaying method for protecting a transformer comprising the steps of: obtaining primary and secondary and voltages and currents of said transformer; calculating flux linkage increments, flux linkage difference increments, the ratio of flux linkage increments of the primary and secondary windings, or the ratio of flux linkage difference increments of the primary and secondary windings from said currents and said voltages; calculating at least one predetermined decision parameter derived from at least one predetermined equation; deciding whether an internal winding fault occurs by comparing said decision parameter to said flux linkage increments, said flux linkage difference increments, said ratio of flux linkage increments of the primary and secondary windings, or said ratio of flux linkage difference increments of the primary and secondary windings; wherein said transformer is a three-phase Y-Δ transformer; and wherein, in the step of calculating at least one predetermined decision parameter derived from at least one predetermined equation, said predetermined equations are ${\frac{{\Delta\lambda}_{A}}{{\Delta\lambda}_{ab}} = \frac{N_{1}}{N_{2}}},\quad {\frac{{\Delta\lambda}_{B}}{{\Delta\lambda}_{bc}} = \frac{N_{1}}{N_{2}}},\quad {\frac{{\Delta\lambda}_{C}}{{\Delta\lambda}_{ca}} = {{\frac{N_{1}}{N_{2}}\quad {or}\quad \frac{\Delta \left( {\lambda_{C} - \lambda_{A}} \right)}{\Delta \left( {\lambda_{ca} - \lambda_{ab}} \right)}} = \frac{N_{1}}{N_{2}}}},{\frac{\Delta \left( {\lambda_{A} - \lambda_{B}} \right)}{\Delta \left( {\lambda_{ab} - \lambda_{bc}} \right)} = \frac{N_{1}}{N_{2}}},\quad {\frac{\Delta \left( {\lambda_{B} - \lambda_{C}} \right)}{\Delta \left( {\lambda_{bc} - \lambda_{ca}} \right)} = \frac{N_{1}}{N_{2}}}$

or modifications of said equations (Δλ_(A), Δλ_(B), Δλ_(C): primary flux linkage increments; Δλ_(ab), Δλ_(bc), Δλ_(ca): secondary flux linkage increments; N₁, N₂: numbers of turns of primary and secondary windings).
 16. The method of claim 15, wherein said decision parameters are: ${{Detector}\quad 1} = {\frac{{\Delta \left( {\lambda_{C} - \lambda_{A}} \right)} - {\frac{N_{1}}{N_{2}}{\Delta \left( {\lambda_{ca} - \lambda_{ab}} \right)}}}{\frac{V_{CArms}}{4.44f}\Delta \quad T} \times 100(\%)}$ ${{Detector}\quad 2} = {\frac{{\Delta \left( {\lambda_{A} - \lambda_{B}} \right)} - {\frac{N_{1}}{N_{2}}{\Delta \left( {\lambda_{ab} - \lambda_{bc}} \right)}}}{\frac{V_{ABrms}}{4.44f}\Delta \quad T} \times 100(\%)}$ ${{Detector}\quad 3} = {\frac{{\Delta \left( {\lambda_{B} - \lambda_{C}} \right)} - {\frac{N_{1}}{N_{2}}{\Delta \left( {\lambda_{bc} - \lambda_{ca}} \right)}}}{\frac{V_{BCrms}}{4.44f}\Delta \quad T} \times 100(\%)}$ (V_(ABrms), V_(BCrms), V_(CArms):primary  rated  voltages).


17. A relaying method for protecting a transformer comprising the steps of: obtaining primary and secondary and voltages and currents of said transformer; calculating flux linkage increments, flux linkage difference increments, the ratio of flux linkage increments of the primary and secondary windings, or the ratio of flux linkage difference increments of the primary and secondary windings from said currents and said voltages; calculating at least one predetermined decision parameter derived from at least one predetermined eguation; deciding whether an internal winding fault occurs by comparing said decision parameter to said flux linkage increments, said flux linkage difference increments, said ratio of flux linkage increments of the primary and secondary windings, or said ratio of flux linkage difference increments of the primary and secondary windings; wherein said transformer is a three-phase Δ-Y transformer; and wherein, in the step of calculating at least one predetermined decision parameter derived from at least one predetermined equation, said predetermined equations are ${\frac{{\Delta\lambda}_{CA}}{{\Delta\lambda}_{c}} = \frac{N_{1}}{N_{2}}},\quad {\frac{{\Delta\lambda}_{AB}}{{\Delta\lambda}_{a}} = \frac{N_{1}}{N_{2}}},\quad {\frac{{\Delta\lambda}_{BC}}{{\Delta\lambda}_{b}} = {{\frac{N_{1}}{N_{2}}\quad {or}\quad \frac{{\Delta\lambda}_{CA} - {\Delta\lambda}_{AB}}{{\Delta\lambda}_{c} - {\Delta\lambda}_{a}}} = \frac{N_{1}}{N_{2}}}},{\frac{{\Delta\lambda}_{AB} - {\Delta\lambda}_{BC}}{{\Delta\lambda}_{a} - {\Delta\lambda}_{b}} = \frac{N_{1}}{N_{2}}},\quad {\frac{{\Delta\lambda}_{BC} - {\Delta\lambda}_{CA}}{{\Delta\lambda}_{b} - {\Delta\lambda}_{c}} = \frac{N_{1}}{N_{2}}}$

or modifications of equations (Δλ_(AB), Δλ_(BC), Δλ_(CA): primary flux linkage increments; Δλ_(a), Δλ_(b), Δλ_(c): secondary flux linkage increments; N₁, N₂: numbers of turns of primary and secondary windings).
 18. The method of claim 17, wherein said decision parameters are: ${{Detector}\quad 1} = {\frac{{\Delta \left( {\lambda_{CA} - \lambda_{AB}} \right)} - {\frac{N_{1}}{N_{2}}{\Delta \left( {\lambda_{c} - \lambda_{a}} \right)}}}{\frac{V_{CArms}}{4.44f}\Delta \quad T} \times 100(\%)}$ ${{Detector}\quad 2} = {\frac{{\Delta \left( {\lambda_{AB} - \lambda_{BC}} \right)} - {\frac{N_{1}}{N_{2}}{\Delta \left( {\lambda_{a} - \lambda_{b}} \right)}}}{\frac{V_{ABrms}}{4.44f}\Delta \quad T} \times 100(\%)}$ ${{Detector}\quad 3} = {\frac{{\Delta \left( {\lambda_{BC} - \lambda_{CA}} \right)} - {\frac{N_{1}}{N_{2}}{\Delta \left( {\lambda_{b} - \lambda_{c}} \right)}}}{\frac{V_{BCrms}}{4.44f}\Delta \quad T} \times 100(\%)}$ (V_(ABrms), V_(BCrms), V_(CArms):primary  rated  voltages).


19. A relaying method for protecting a transformer comprising the steps of: obtaining primary and secondary and voltages and currents of said transformer; calculating flux linkage increments, flux linkage difference increments, the ratio of flux linkage increments of the primary and secondary windings, or the ratio of flux linkage difference increments of the primary and secondary windings from said currents and said voltages; calculating at least one predetermined decision parameter derived from at least one predetermined equation; deciding whether an internal winding fault occurs by comparing said decision parameter to said flux linkage increments, said flux linkage difference increments, said ratio of flux linkage increments of the primary and secondary windings, or said ratio of flux linkage difference increments of the primary and secondary windings; wherein said transformer is a three-phase Δ-Δ transformer; and wherein, in the step of calculating at least one predetermined decision parameter derived from at least one predetermined equation, said predetermined equations are ${\frac{{\Delta\lambda}_{CA}}{{\Delta\lambda}_{ca}} = \frac{N_{1}}{N_{2}}},\quad {\frac{{\Delta\lambda}_{AB}}{{\Delta\lambda}_{ab}} = \frac{N_{1}}{N_{2}}},\quad {\frac{{\Delta\lambda}_{BC}}{{\Delta\lambda}_{bc}} = {{\frac{N_{1}}{N_{2}}\quad {or}\quad \frac{{\Delta\lambda}_{CA} - {\Delta\lambda}_{AB}}{{\Delta\lambda}_{ca} - {\Delta\lambda}_{ab}}} = \frac{N_{1}}{N_{2}}}},{\frac{{\Delta\lambda}_{AB} - {\Delta\lambda}_{BC}}{{\Delta\lambda}_{ab} - {\Delta\lambda}_{bc}} = \frac{N_{1}}{N_{2}}},\quad {\frac{{\Delta\lambda}_{BC} - {\Delta\lambda}_{CA}}{{\Delta\lambda}_{bc} - {\Delta\lambda}_{ca}} = \frac{N_{1}}{N_{2}}}$

or modifications of said equations (Δλ_(AB), Δλ_(BC), Δλ_(CA): primary flux linkage increments; Δλ_(ab) , Δλ_(bc), Δλ_(ca): secondary flux linkage increments; N₁, N₂: numbers of turns of primary and secondary winding).
 20. The method of claim 19, wherein said decision parameters are: ${{Detector}\quad 1} = {\frac{{\Delta \left( {\lambda_{CA} - \lambda_{AB}} \right)} - {\frac{N_{1}}{N_{2}}{\Delta \left( {\lambda_{ca} - \lambda_{ab}} \right)}}}{\frac{V_{CArms}}{4.44f}\Delta \quad T} \times 100(\%)}$ ${{Detector}\quad 2} = {\frac{{\Delta \left( {\lambda_{AB} - \lambda_{BC}} \right)} - {\frac{N_{1}}{N_{2}}{\Delta \left( {\lambda_{ab} - \lambda_{bc}} \right)}}}{\frac{V_{ABrms}}{4.44f}\Delta \quad T} \times 100(\%)}$ ${{Detector}\quad 3} = {\frac{{\Delta \left( {\lambda_{BC} - \lambda_{CA}} \right)} - {\frac{N_{1}}{N_{2}}{\Delta \left( {\lambda_{bc} - \lambda_{ca}} \right)}}}{\frac{V_{BCrms}}{4.44f}\Delta \quad T} \times 100(\%)}$ (V_(ABrms), V_(BCrms), V_(CArms):primary  rated  voltages). 